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GENERATION OF A DESIRED WAVEFRONT WITH A PLURALITY OF PHASE 
CONTRAST FILTERS 

FIELD OF THE INVENTION 

The present invention relates to a method and a system for synthesizing a prascribed 
intensity pattem based on generalized phase contrast imaging. 
BACKGROUND OF THE INVENTION 

It is well known to form an Image by phase contrast Imaging methods in which phase 
modulation of light is converted into Intensity modulation. As opposed to Intensity 
modulation, phase modulation does not Involve loss of energy. 

A generalized phase contrast imaging method and system is disclosed in International 
Patent Appl.cat,on Publication No. WO 96/34207 for synthesizing a prescribed Intensity 

the phase shrft 0 ,s less than 1 radian. An improved method is provided without this 
assumption and based on an imaging operation with a simple one-to-one mapping between 
rasou^on elements or pixels of a spatial phase modulator and resolution element! of the 
generated intensity pattern. 

di^osed phase conh^s, ,™ging method of syn«,es*9 3 p,^bed intensity pattern 
i(x,y),compnses the steps of: 

• dividing the intensity pattem l(x'y) Into pixels In accordance with the disposition of 
resolution elements (x.y) of a spatial phase mask having 

a plurality of Individual resolution elements (x.y). each resolution element (x y) 
modulating the phase of electromagnetic radiation Incident upon it with a ' 
predetennlned phasor value e'*^^, 

• radiating electromagnetic radiation towards the spatial phase mask. 

• Fourier or Fresnel transforming the modulated electromagnetic radiation. 

. Phase shifting with a spatial phase filter (SPF) In a region of spatial frequencies 

comprising DC In the Fourier or Fresnel plane, the modulated electromagnetic radiation 
by a predetemiined phase shift value e in relation to the remaining part of the 
electromagnetic radiation, and 

. fonming the Inlensl.), pattern by Fourier or Fresnel t^nsformlng, respectively, the phase 
shrfted Founeror Fresnel transfonned modulated elect™n«g„et,o radiation, whereby 
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each resolution element (x,y) of the phase mask Is imaged on a corresponding 
resolution element (x'.y') of the Image, 

. calculating the phasor values e^-^^ of the phase mask and the phase shift value 9 in 
accordance with 

for selected phase shift values B. a being the average of the phasors of the 
resolution elements of the phase mask, 

• selecting, for each resolution element, one of two phasor values which represent a 
particular grey level, and 

. supplying the selected phasor values e'«^v) to the resolution elements (x,y) of the spatial 
phase mask. 

In one embodiment disclosed In WO 96/34207, the spatial phase mask is positioned at the 
front focal plane of a lens while the spatial phase filter Is positioned In the back focal plane of 
the lens, whereby a first electromagnetic field at the spatial phase mask is Fourier 
transformed by the lens into a second electromagnetic field at the spatial phase filter. Thus 
specific spatial frequencies of the first electromagnetic field will be transmitted through the ' 
spatial phase filter at specific positions of the phase filter. For Instance, the energy of the 
electromagnetic radiation at zero frequency (DC) Is modified by the phase filter and 
transformed onto the Intersecting point of the Fourier plane and the optical axis of the lens 
also denoted the zero-order diffraction region. 

Typically, the spatial phase filter is adapted to phase shift the DC-part of the electromagnetic 
radiation and to leave the remaining part of the electromagnetic radiation unchanged 
Alternatively, the DC-part of the electromagnetic radiation remains unchanged and the 
remaining part of the electromagnetic radiation is phase shifted. This alternative is preferred 
when the energy level of the DC-part of the electromagnetic radiation is so high that the 
phase shifting part of the phase filter will be destroyed by It. For example in laser cutting the 
DC level of the laser beam can be so high that a phase shifting dot positioned at the 
intersecting point of the DC part of the laser beam at the phase filter would deteriorate It is 
also possible to block the electromagnetic radiation (no ti^nsmittance) in the zero-order 
diffraction region, however, the DC energy of the radiation Is then lost. 

In US 2003/0030902 a microscope system is disclosed for obtaining images of optimum 
image quality comprising 1) an illuminating device which emits light from a light source to an 
object, and generates a luminous fiux Including infomiation of the object. 2) an Illuminating 
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Hght modulating device which modulates at least one of wavelength, phase, intensity 
polarization, and coherency of the light emitted to the object by the illuminating device 3) an 
objective lens and an imaging lens which focus the luminous flux including the information of 
the object to form the image of the object. 4) a pupil modulating device which is disposed 
near a pupil plane of the objective lens, and modulates at least one of phase, intensity and 
direction of polarization of the luminous flux including the infomiation of the object 5) an 
Image pickup device which is disposed on a plane on which the image of the object is 
formed by the objective lens and the imaging lens, and picks up the image of the object 6) 
an image display device which displays the image of the object picked up by the image 
pickup device. 7) an image analysis device which analyses the image of the object picked up 
by the image pickup device, and 8) a parameter decision device which adjusts the 
modulation amounts of the illuminating light modulating device and the pupil modulating 
device by using the image infomiatlon of the object analysed by the image analysis device. 
The pupil modulating device is a phase filter and in one of the disclosed embodiments (Fig. 
7) the phase shift is equal to | and the phase shifting area of the phase filter is annular. In 
another of the enclosed embodiments (Fig. 9) the phase shift is equal to n and the phase 
filter IS divided into a plurality of rows, every second row having phasor value e^ and being 
interlaced with the remaining rows having the phasor value e'° 

US 2003/0030902 does not disclose a method of synthesizing an intensity pattern by 
designing a phase modifying element with phasor values e^^ calculated from Intensity 
values Of the desired intensity pattern. In US 2003/0030902. the phase modifying element is 
a sample that is obtained for studying It in a microscope. The sample is not designed with 
the purpose of generating a certain Intensity pattern. 
SUMMARY OF THE INVENTION 

The previously mentioned known method of synthesizing a specific Intensity pattern 
operates on a plane incoming electromagnetic field. 

It is an Object of the present invention to provide a method and a system for synthesizing an 
intensity pattern operating on a more complex electromagnetic field than the known method 
for further flexibility, e.g.. to be able to utilize light sources that are not based on a single 
collimated point source, such as laser arrays. VCSEL-an-ays. LED-arrays etc. 
Further, the spatial phase filter according to the known method processes an incident 
electromagnetic field In the very small diffraction region only. e.g. the zero-order diffraction 
region, leaving the remaining part of the filter un-utllized. Distributing the spatial processing 
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of an electromagnetic field across a larger area of the spatial phase filter would make the 
filter capable of handling higher energy levels. Further, it should be possible to match the 
phase shifting regions of the spatial phase filter to a source for emission of electromagnetic 
radiation that cannot be considered to constitute a point source, because it emits radiation 
from a large area. 

According to a first aspect of the present invention, the above and other objects are fulfilled 
by a phase contrast system for synthesizing an Intensity pattern l(x'. /). comprising 

• a source of electromagnetic fields for emission of at least two substantially plane 
electromagnetic fields with different axes of propagation, 

• a phase modifying element for phase modulation of the electromagnetic fields by 
Phaser values e^^^ and positioned so that the at least two electromagnetic fields are 
incident upon it at different respective angles of approach. 

• first Fourier or Fresnel optics for Fourier or Fresnel transforming the.phase modulated 
electromagnetic fields positioned in the propagation paths of the at least two phase 
modulated fields, 

• a spatial phase filter with at least two phase shifting regions positioned at respective 
zeroKDrder diffraction regions of the at least two respective phase modulated 
electromagnetic fields for individually phase shifting the at least two respective Fourier 
or Fresnel transfomied electromagnetic fields by predetermined respective phase shift 
values e„ in relation to the remaining part of the at least two respective transfomied 
electromagnetic fields, and 

» second Fourier or Fresnel optics for fom^ing the intensity pattern I(x'. /) by Fourier or 
Fresnel transfom^ing the at least two respective phase shifted Fourier or Fresnel 
transformed electromagnetic fields, 

the Phaser values e'«^» ef the phase modifying element and me phase shift values e„ 

substantially fulfilling that 

^(^^;'0^Z5(«K|exp(4(x^;;0) + ^Ja|K^-•exp(z•e„)-l)f 
for selected phase shift values 9„, wherein 

A is an optional amplitude modulation of the spatial phase filter outside the zero-order 
diffraction regions, 
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Bn is an optional amplitude modulation of the spatial phase filter in the respective n'th 
zero-order diffraction region, 

a = Hexp(/j^) Is the average of the phasors e^'^ of the resolution elements of the 
phase modifying element, and 

^ = and 

S{n) is the Intensity of the n'th electromagnetic field, and 

= 1 - (l ■227tTj„ ) , wherein 
Jo is the zero-order Bessel function and 

% relates the radius R,„ of the n'th zero-order filtering region to the radius R, of the 
main-lobe of the Airy function of the Input aperture, r]„=R^/Jt^. 

According to a second aspect of the present invention, the above and other objects are 
fulfilled by a method for synthesizing an Intensity pattern l(x'. /). comprising the steps of 
. dividing the intensity pattern l(x'./) into pixels In accordance with the disposition of 
resolution elements (x.y) of a phase modifying element having 

a plurality of individual resolution elements (x.y), each resolution element (x y) 
modulating the phase of electromagnetic radiation incident upon it with a 
predetermined phasor value e'*^'*''^, 

. radiating at least two substantially plane electromagnetic fields with different axes of 
propagation towards the phase modifying element so that the at least two 
electromagnetic fields are incident upon it at different respective angles of approach. 

. Fourier or Fresnel transfomiing the phase modulated electromagnetic fields, 

. Phase Shifting In at least two phase shifting regions positioned at respective zero-order 
d.ffract.on regions of the at least two respective phase modulated electromagnetic fields 
for ,nd,v.dually phase shifting the at least two respective Fourier or Fresnel transformed 
electromagnetic fields by predetemilned respective phase shift values 9. in relation to 
the remaining part of the at least two respective transformed electromagnetic fields, and 

• fomiing the intensify pattern I(x-, y) by Fourier or Fresnel transforming the at least two 
respective phase shifted Fourier or Fresnel transfomied electromagnetic fields. 
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• calculating the phasor values e'^^ v) of the phase modifying element and the phase shift 
values e„ substantially in accordance with 

^(^^J'0 = S•S(«Kiexp(4(^^;'')) + ^„|a|(B„^-•exp(^e„)-l)p 
for selected phase shift values 9„. wherein 

A is an optional amplitude modulation of the spatial phase filter outside the zero- 
order diffraction regions, 

Bn is an optional amplitude modulation of the spatial phase filter in the respective n'th 
zero-order diffraction region, 

a = \a\ exp {i^ ) is the average of the phasors e^^^ of the resolution elements of the 
phase modifying element, and 

^ = and 

Sin) is the intensity of the n'th electromagnetic field, and 

■^n = 1 - ^0 (1 ■22jrTj„ ) , wherein 

Jo is the zero-order Bessel function and 

Vn relates the radius of the n'th zero-order filtering region to the radius of the 
main-lobe of the Airy function of the input aperture, ti„ = i?i„ /ilj . 

. selecting, for each resolution element, one of two phasor values which represent a 
particular grey level, and 

. supplying the selected phasor values e'^'^v) to the respective resolution elements (x. y) 
of the phase modifying element. 

The mathematical expressions will be further explained below. 

The axis of propagation of a plane electromagnetic field is perpendicular to the electric and 
magnetic fields. 

In an embodiment of the present Invention, the phase modifying element has an input for 
reception of signals for addressing the resolution elements (x. y) and for adjusting the 
phasor values e^-^ of the respective addressed resolution elements (x. y). 
The system may further comprise a controller with a first output that is connected to the input 
of phase modifying element, and a second output that is connected to the spatial phase filter 



wo 2004/113993 

PCT/DK2004/000452 

-7- 

and being adapted for adjusting phasor values e'^^ v) of the phase modifying element and 
phase shift values 9„ of the spatial phase filter. Preferably, the at least two substantially 
plane electromagnetic fields are non-interfering i.e. the fields are mutually Incoherent and/or 
orthogonally polarized and/or have different wavelengths, etc. 

It should be noted that, in each resolution element of the phase modifying element, one of 
two phasor values which represent a particular grey level of the Intensity pattem l(x'.y) may 
be selected. 

In an embodiment of the present invention, the spatial phase filter substantially does not 
attenuate the electromagnetic fields incident upon It outside the phase shifting regions. 1. e. 
A is equal to one or approximately equal to one. 

In an embodiment of the present Invention, the spatial phase filter substantially does not 
attenuate the electromagnetic fields incident upon it inside the phase shifting regions, i.e. B„ 
is equal to one or approximately equal to one. 

It is also preferred that the phase shift value substantially fulfils the equation 

1 



2|sin^ 
2 



for a loss less filter with ^ = 1 and 5„ = 1 .in a preferred embodiment of the present 
Invention, the phase shift e„ is equal to nor approximately equal to n. Accordingly the 
previous equation leads to K„a = X and the phase values H^.yU the phase modifying 
element may be calculated in accordance with 

A 

K^A-' l\sm(<l>{x,y))dxdy = 0 
where /I Is the illuminated area of the phase modifying element. 

The electromagnetic field or radiation may be of any frequency range of the electromagnetic 
spectrum, i.e. the gamma fi-equency range, the ultraviolet range, the visible range the 
infrared range, the far Infrared range, the X-ray range, the microwave range, the HF (high 
frequency) range, etc. The present Invention is also applicable to particle radiation, such as 
electron radiation, neutron radiation, etc. 

Preferably, the electromagnetic fields are monochromatic or quasi-monochromatic so that 
the energy of the electromagnetic fields is concentrated in a narrow frequency bandwidth. 
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Since the intensity pattern is reconstructed by. for each of the Incoming plane fields. 
Interference of two electromagnetic fields generated by different phase shifting of different 
parts of the Incoming field, it is required that the frequency range of the emitted 
electromagnetic field is sufficiently narrow to ensure that the two electromagnetic fields are 
coherent so that their superposition generates the desired Intensity pattern. If the frequency 
range is too broad, the two fields will be incoherent and the phase Information will be lost as 
superposition of non-coherent fields results In a summation of the intensities of the two 
fields. It is required that the difference between individual delays of electromagnetic fields to 
be superposltloned Is less than the wavelength of the fields. This Is a relaxed requirement 
that allows the electromagnetic fields to be relatively broad-banded. For example in the 
visible range, a Xe-lamp or a Hg-lamp can be used as a light source in a system according 
to the present Invention with the advantage compared to a laser light source that speckle 
noise is reduced. The requirements of the spatial coherence of the electromagnetic fields 
depend upon the space bandwidth product of the con-esponding system and how close the 
required system perfomiance Is to the theoretically obtainable perfomiance of the system. 
Preferably, the electromagnetic radiation is generated by a coherent source of 
electromagnetic radiation, such as a laser, a semi-conductor laser, a strained multi-quantum 
well laser, a vertical cavity surface emitting laser (VCSEL). a maser. a phase-locked laser 
diode array, a light emitting diode, a pulsed laser, such as a sub-picosecond laser, etc. or an 
array of such sources. However, as already mentioned, a high-pressure arc lamp, such as 
an Hg lamp, a Xe lamp. etc. may also be used and even an incandescent lamp may be used 
as a source of electromagnetic radiation. 

The phase modifying element changes the phase of an electromagnetic field incident upon 
it. Optionally, it may also change the amplitude of an electromagnetic field Incident upon it. 
The phase modifying element may transmit or reflect the incident electromagnetic field. The 
phase modifying element is divided into a number of resolution elements, each of which 
modulates the incident electromagnetic field by changing its phase by a specific 
predetermined value. The predetermined values are assigned to each resolution element in 
different ways depending upon the technology applied in the component. For example in 
spatial light modulators, each resolution element may be addressed either optically or 
electrically. The electrical addressing technique resembles the addressing technique of 
solid-state memories in that each resolution element can be addressed through electronic 
circuitry to receive a control signal corresponding to the phase change to be generated by 
the addressed resolution element. The optical addressing technique addresses each 
resolution element by pointing a light beam on it, the Intensity of the light beam 
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corresponding to the phase change to be generated by the resolution element illuminated by 
the light beam. 

Spatial phase modulation may be realized utilizing fixed phase mask, a liquid crystal device 
based on liquid crystal display technology, a MEMS (micro electro-mechanical system), a 
MOEMS (micro opto-electro-mechanical system), such as a dynamic mirror device, a digital 
micro-mirror an-ay. a defomiable mirror device, etc, a membrane spatial light modulator, a 
laser diode array (integrated light source and phase modulator), smart pixel arrays, etc' 
Seiko-Epson produces a transmitting liquid crystal SLM (LC-SLM) having a high resolution 
matrix of transparent liquid crystal elements wherein the relative permitbvity of each element 
can be electrically modulated in order to vary the refiractive Index and thereby the optical 
path length of the element. 

Meadowlark produces a parallel-aligned liquid crystal (PAL-SLM) with a high fill factor, but 
this device has a very low resolution In that it contains only 137 phase-modulating elements. 

Hamamatsu Photonics produces a dynamically controllable PAL-SLM with VGA or XGA 
resolution. 

Texas Instruments produces a Digital Min-or Device (DMD) having an array of mirrors each 
of which can be tilted between two positions. 

In an embodiment of the invention, the phasor values are contained in an object that Is not 
controlled by the system, e.g. a specimen In a slide, in stead they are detected in the 
system. 

The spatial phase filter Is typically a fixed phase mask, such as an optically flat glass plate 
coated with a dielectric layer in the region wherein the modulated electromagnetic field is 
phase shifted in relation to the remaining part of the electromagnetic field. However, the 
spatial phase modulators mentioned in the previous section may also be used for spatial 
phase filters. In addition, non-linear materials providing for self-phase modulation such as 
Ken--type materials can also be used for introducing phase shift 

An imaging system maps the phase modulating resolution elements of the phase modifying 
element on the target surface of the reconstmcted intensity pattem. It may comprise a 4f- 
lens configuration (two Fourier transfomiing lenses utilizing transmission of light or one 
Fourier transfomiing lens utilizing reflection of light) or a single imaging lens. However, any 
optical imaging system providing a filtering plane for the spatial phase filter may be applied 
in a phase contrast imaging system. 
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genemted bysupen=osi«on CeleCromagnetlc fields ,„ me image plane ofthe Imaging 

The Phase modifying element changes ,he phase values of an eleotromagnete field 
.noden. upon ,. and Imaging system directs the electn^magnetlc «e,d with chafed 
Phases ^fleeted from ort:ansmltted through the phase modifying element towards the 
pa*al phase Alter. The phase Alter phase shifts a part of «,e electromagnetic field and the 
imaging system ,s adapted to supenmpose in the Image plane the phase shifted part of the 

According to a preferred embodiment of the Invention, the phase modifying element is 
Portioned a, the front foca, plane of a lens while the spaAai phase Alter is posiAoned in Ate 
TIT T- . " ^'-^"'-aneAc Aeld a, the phase modifying 

e^em ,s Founer transformed by the lens into a second elecAomagneAc fleid a. Jphase 

mrough the spatial phase Alter at speclAc poslAons of *,e phase fliter. For Instance the 

ZZ "^""^ *e Phase 

fliter at the .ntersecAng point of the Founer plane and the optical axis of A,e lens also 

denoted the zero-order diffraction region. 

in 'Ztr'^T °' °' ^ Of source ,s faalitated 

matched to A,e geometry of Are source. For example, if a linear arrw of VCSE1.S fomrs L 
source the phase shifting regions of A« spaAal phase Alter fdm, a conespondlng linear 
ar^y of phase shWng regions, each of A« regions being poslAoned at Are zere^rder 
drffracAon region of a respecOve VCSEL In the VCSEL array. FurAner, A,e shape of each 
phase shiftrng region may match Are shape of A» zer^rxJer diffracAon region of Are 
respective VCSEL. 

Ukewlse, a phase fliter may match a source v*, a speclAc geometrical shape wlA, a 
conAnuous phase shifting region covering an area of Are phase fliter mat conesponds to Are 
zeroorderdlffracAon region of the source. sponostottre 

Thus. A,e energy of me electmmagnetlc fields of the system may be dlsWbuted over a large 
area compared to me ar«a of a zeroK>rder diffraction region of a single plane 
electromagneAe flekJ of a known phase contrast imaging system. 
Thus the phase shifting regions of the spaBal phase fliter may fbm, a rectangular array a 
aroular anay, a linear array, two linear crossing anays. a conAnuous region, a ring, etl' 
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The at least two substantially plane electromagnetic fields with different axes of propagation 
may be generated in a time multiplexed manner, e.g. by a scanning mirror or prism, 
deflecting or reflecting a beam of electromagnetic field in different directions of propagation. 
The capability of handling high energy levels of electromagnetic fields of the present 
invention may be utilized for provision of a laser cutter comprising a system according to the 
present invention. 

Further, the capability of handling high energy levels in combination with the capability of 
generating crossing light beams may be utilized for provision of an optical tweezer or an 
anay of optical tweezers according to the present Invention. 

In an embodiment of the present invention, wherein the apertures of the system is of 
insignificant importance to the operation of the system and calculation of the phasor values. 
I.e. Kn Is equal to one or approximately equal to one. 

For a more detailed understanding of the invention, the Zemike approximation is reviewed 
below, followed by a generalization where the above-mentioned mathematical expressions 
are derived for a single (n=1) on-axis centred phase contrast filtering implementation. 
The Zemike phase contrast method allows for the visualization of phase perturbations by the 
use of a Fourier plane phase shifting filter. The Dutch physicist Fritz Zemike received the 
Nobel Prize in 1953 for Inventing this method, which led to a break-through in medicine and 
biology by making essentially transparent cell or bacteria samples clearly visible under a 
microscope. Its successful operation, however, requires that the spatial phase distribution. 
^x,y) , at the input Is limited to a "small-scale" phase approximation where the largest 
phase is typically taken to be significantly less than ;r / 3 . If the phase distribution at the Input 
is thus restricted, then a Taylor expansion to first order is sufficient for the mathematical 
treatment so that the Input wavefront can be written as 

exp (i^ (x, j;)) « 1 + z<|) (x, y) (1 ) 

The light con-esponding to the two terms In this "small-scale" phase approximation can be 
separated spatially by use of a single lens where the phase distribution is located in the front 
focal plane and the con-esponding spatial Fourier transformation is generated in the back 
focal plane of the lens. In this first order approximation the constant term represents the 
amplitude of on-axis light focused by the lens In the back focal plane and the second 
spatially varying term represents the off-axis light. Zemike realized that a small phase 
shifting quarter wave plate acting on the focused light makes it possible to obtain an 
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approximately linear visualization of small phase structures by generating interference 
between the two-phase quadrature terms in Eq. (1): 

/(*'.>'')«l+2<|)(*',:»/') (2) 

it Should be noted that a three-quarter waveplate works equally well to produce contrast but 
the plus sign in Eq. (2) is negated leading to so-called negative phase contrast. A substantial 
improvement in the visibility of the Zemike phase contrast visualization In Eq. (2) requires 
strong damping of the focused light in addition to the phase shift required to generate the 
contrast. 

In the general case, where we are not limited to a small-scale Input phase perturbation we 
cannot assume that a series expansion to first order as in the Zernike approximation Is a 
sufficient representation of a given phase perturbation. Higher order terms In the expansion 
need to be taken into account, so the expansion takes the form: 

exp(4(^.J')) = l + '<|.(^.;')-i<l>M^.J')-i'V(*.>') + ^r(;c.;;)+... (3) 

However, here the spatially varying terms can not be considered as separate from the 
supposedly focused light represented by the first term in this Taylor series expansion as is 
implied by the Zemike approach, and all of these spatially varying temis contribute to the 
intensity of the on-axis focused light. For a significant modulation in the input phase this 
contnbutlon of the spatially varying temis can result In a significant modulation of the focal 
spot amplitude In the back focal plane of the lens. These temis can in fact result in either 
constructive or destructive Interference with the on-axis light, although the net result will be 
an attenuation of the focused light amplitude, which only has a maximum value for a perfect 
unperturbed plane wave at the Input. 

For phase objects breaking the Zemike approximation we must, therefore, find an alternative 
mathematical approach to that of the Taylor expansion given in Eq. (3). We have chosen a 
Fourier analysis as a more suitable technique for completely separating the on-axis and 
higher spatial frequency components. This gives the following fonn for exp{i^x,y)), where 



M'H^'y)) = f JWl jjMM^'y))dxdy + -higher ftequency terms" 



(4) 



In this Fourier decomposition the first tem, Is a complex valued constant linked to the on-axis 
focused light from a phase object defined within the spatial region. D. and the second temi 
descnbes light scattered by spatially varying structures In the phase object. Comparing Eq 
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(3) and Eq. (4) it is apparent that the first term of Eq. (3) is a poor approximation to the first 
term of Eq. (4) when operating beyond the Zemil<e small-scale phase regime. 
An important issue to consider when analysing the effect of spatial filtering on the light 
diffracted by phase perturbations is the definition of what spatially constitutes focused and 
scattered light. In the previous description of Zemike phase contrast It was assumed that the 
focused light is spatially confined to a somewhat unphysical delta function. As we know any 
aperture truncation such as occurs within a practical optical system will lead to a 
corresponding spatial broadening of the focused light It Is therefore essential that we define 
the terms "focused iighf and "scattered light" explicitly for such a system. In this context It Is 
necessary to look more carefully at the sequence of apertures confining the light wave 
propagation through a typical optical set-up. 

A commonly applied architecture that provides an efficient platform for spatial filtering Is 
Illustrated in Fig. 1 and is based on the so-called 4-f configuration. An output interferogram 
of an unknown phase object or phase disturbance is obtained by applying a truncated on- 
axis filtering operation In the spatial frequency domain between two Fourier transforming 
lenses (LI and 12). The first lens performs a spatial Fourier transform so that directly 
propagated light Is focused Into the on-axis filtering region whereas spatially varying object 
mformation generates light scattered to locations outside this central region. We can 
describe a general Fourier filter in which different phase shifts and amplitude damping 
factors are applied to the "focused" and "scattered" light. In Fig 1. we show a circularly 
symmetric Fourier filter described by the amplitude transmission factors A and 8 for the 
"scattered" and "focused" light respectively and by the relative phase shift ft These filter 
parameters can be chosen to replicate any one of a large number of commonly used filter 
types 0.e. phase contrast. dart< central ground, point diffraction and field-absorption filtering) 
By applying a given Fourier filter and a second Fourier lens, we obtain an interference 
pattern in the observation plane. The focused on-axis light acts as the synthetic reference 
wave (SRW) in the CPI system, this interferes with the scattered light to generate the output 
interference pattern. In the following section we discuss the Importance of the SRW and 
show how it influences, among other things, the choice of the Fourier filter parameters. 
Having described the generic optical system that makes up the CPI. we turn to a detailed 
analytical treatment of the Important elements in this system. Assuming a circular input 
aperture with radius. Ar . truncating the phase disturbance modulated onto a colllmated unit 
amplitude, monochromatic field of wavelength. ;i . we can describe the incoming light ' 
amplitude a(x,y) by. 
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« (jcj')=circ(//Ar ) exp (i^ (x, y)) (5) 
at the entrance plane of the optical system shown in Fig. 1 using the definition that the drc- 
functlon is unity within the region, r = V^T/ ^ Ar . and zero elsewhere. 
Similarly, we assume a circular on-axis centred spatial filter of the form: 

^(/,./,) = ^[l + (5^-'exp(,e)-l)circ(/,/A/,)] (g) 

where B e [O; l] is the chosen filter transmittance of the focused light. 0 e [0;2;r] is the 
applied phase shift to the focused light and A 6[0;l] is a filter parameter describing field 
transmittance for off-axis scattered light as indicated in Fig. 1 . The spatial frequency 
coordinates are related to spatial coordinates in the filter plane such that: 
(f.>f,) = i¥r\x,,y^) and =4fFTf[. 

Perfomiing an optical Fourier transform of the input field from Eq. (5) followed by a 
multiplication with the filter parameters in Eq. (6) and a second optical Fourier transform 
(con-esponding to an inverse Fourier transform with inverted coordinates) we obtain an 
expression for the intensity /(x',y) describing the interferogram at the observation plane of 
the set-up: 

K^'.y) = ^'|exp(/^(j:',:,;'))circ(ryAr) + ^BA'' exp(/^) - l)g(r')\' (7) 

where ^(r') is the synthetic reference wave (SRW) and the terms a and -^(x'.y) are given 
by: 

a=(7t(Ar)') ' JJ exp(/,|,(x,j^))iaf^ = |a|exp(«t.-) 

(8) 

it should be noted that to achieve a tractable analytic expression in Eq. (7) it has been 
assumed that the spatial frequency content of the phase object is sufficiently described by 

the term, a . within the on-axis centred filtering region characterized by the spatial frequency 
range A/, . 

The generally complex valued and object dependent term, a , corresponding to the 
amplitude of the focused light plays a significant role in the expression for the interference 
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pattern described by Eq. (7). Referring to the discussion in the introduction, we are now able 
to confirm that the frequent assunnption, that the amplitude of the focused light is 
approximately equal to the first term of the Taylor expansion in Eq. (1 ). can generally result 
in misleading interpretations of the interferograms generated at the CPI output. 

Of similar importance in the analysis of Eq. (7) is the temi g{r') describing the spatial profile 
of the SRW. diffracted from the aperture formed by the on-axis centred filtering region. It is 
the interference between this SRW term, carrying the information about the filtering 
parameters, and the Imaged phase object that generates the output Interferogram. Thus it is 
important to obtain an accurate description for the SRW and thereby an accurate derivation 
for Eq. (7). The zero-order Hankel transform followed by a series expansion in the spatial 
dimension, r'. will be used to describe the SRW. This is a relatively simple approach, which 
to the knowledge of the author has not previously been applied to this problem. 
For a circular input aperture with radius. Ar . we can describe the radius of the 
corresponding central phase shifting region of the Fourier filter (characterized by the 
parameters 8 and 6f) in temis of a radial spatial frequency range A/; . We can thus obtain 
the following expression for the SRW by use of the zeroorder Hankel transform: 

^(r') = 2«Arfv,(27rA^,)/,(2nr'/,)4r, (g) 

In order to simplify the analysis, we introduce a term 7. which explicitly relates the radius of 
the central filtering region, i?, , to the radius of the main-lobe of the Airy function. . 
resulting from the Fourier transform of the circular Input aperture alone. We can thus 
express 7 in terms of Ar and A/; such that: 

•n = i?,/^=(0.6l)-VAr, (10) 

where the factor of 0.61 arises from the radial distance to the first zero crossing of the Airy 
function conesponding to half of the Airy mainlobe factor, of 1 .22. If we make this 
substitution in Eq.(9) and then perfomi a series expansion in r' . we obtain the following 
expression for the SRW: 



si'-') = 1-^0 (l.227rn)-[(0.6l7:7i)' J, (l.227CTi)](r 7 Ar)' + 
{[(0.6l7ni)V4][2y3(l.227iTi)-0.61n7iy,(l.22jiTi)]|(ryAr)\. 



(11) 
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In this expansion, the SRW is expressed in radial coordinates normalised to the radius of the 
.maged Input aperture. This can easily be scaled to allow for a magnification within the 
.maging system, though for the remainder of our analysis a direct imaging operation is 
assumed. From Eq.(1 1 ) it is apparent that the SRW will change as a function of the radius of 
the central filtering region. Additionally, It is clear that the SRW profile is not necessarily flat 
over the system output aperture. This is an important, yet often neglected, factor in 
detennining the performance of a CPI. 

Depending on the accuracy needed for the description of the Interferograms one can choose 
to include a number of spatial higher order terms from the expansion In Eq. (1 1 ). The 
.nfluence of the higher order terms has the largest impact along the boundaries of the 
.maged aperture. For rj -values smaller than 0.627 and when operating within the central 
region of the image plane, spatial higher order temis are of much less significance and we 
can safely approximate the synthetic reference wave with the first and space invariant term 
InEq. (11): 

g (r • e central region) « 1 - y„ (l .22%f\) (12) 
SO that we can simplify Eq. (7) to give: 

I{x\y) = A^\^{i^x\y)yK^[BA-' exp(zG)-l)|' (13) 

where K=\-J,{l227C7j). The influence of the finite on-axis filtering radius on the focused 
light is thus effectively included as an extra "filtering parameter so that the four-parameter 
filter set (A,B,9,K{rj)) together with the complex object dependent term, a . effectively 
defines the type of filtering scheme we are applying. 

Having detemilned a suitable operating range for the CPI in terms of the production of a 
good SRW, we must now examine the role that the remaining filter parameters play in the 



From Eq. (13) we see that the filter parameters {A^B,e) can be combined to form a single 
complex valued term. C, the combined filter term, such that: 

C = |C|exp(iV^) = BA-' exp(/e)-l (14) 
therefore. Eq. (13) can be simplified to give: 
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= ^|l+2\C\cos{vrc)+\cf 
sin6? = (5^-')"'|C|sin(t/c) 



(16) 



Since It is a complex variable, the combined filter temi C. which effectively describes the 
complex filter space, can be considered to consist of a vector of phase and length |c| 
as shown In Eq.(14). Thus In order to obtain an ovemew of the operating space covered by 
all the possible combinations of three independent filter parameters {A,B,0) we can now 
instead choose to consider a given filter In temis of the two combined parameters and 
|C| . However, refemng to Eq.(15). It can be seen that the filter parameter. A, also appears 
independently of the combined filter term. C. Fortunately, this Issue can be resolved by 
considering that the tem, BA'^ from Eq. (14) must be constrained In the following way: 

SA-'<1 => A = l B=\C + 1\ 

BA~^=1 => A = \ R-i 

^' (17) 
BA-'>1 => B = l, ^ = jc + lj-' 

These constraints arise from the adoption of a maximum irradiance criterion minimising 
unnecessary absorption of light In the Fourier filter, which reduces both in-adiance and the 
signal to noise ratio In the CPI output. 

In the prBvlous sections we derived expressions and graphical analysis methods relating the 
spatial average value of a given phase disturbance to obtain peak irradiance and optimal 
visibility in combination with high accuracy in systems with unknown wavefront phase 
disturt^ances. We saw that if a CPI Is applied to wavefront sensing or the visualisation of 
unknown phase objects the Generalised Phase Contrast (GPC) method specifies the filter 
phase and aperture size parameters for achieving optimal performance In extracting and 
displaying the phase information earned by the incoming wavefront. On the other hand In 
cases Where we have control over the incoming wavefront or phase modulation the GPC 
method provides extra means of optimisation by encoding the phase distribution Itself in 
addition to modifying the filter parameters. The two main scenarios: A) synthesizing the 
spatial phase for intensity display or B) measuring the spatial phase with high accuracy 
strongly Influences which of the parameters in the analysis that should be kept fixed and 
Which could be Changed or adapted. The first approach is particularly useful when the filter 
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parameters have a restricted dynamic range or are fixed. The rigorous derivation of the 
equations for choosing these parameters wlil be derived in this section. 

When synthesizing an input phase distribution for optimal visibility of an output Intensity 
pattem the situation Is more relaxed than the situation Involving accurate interferometric 
measurements of unknown phase disturbances. The parameter 77 can therefore In most 
cases be chosen to completely encompass the zero-order light with the result that the term. 
K , tends to unity as the Bessel function tends to zero In Eq. (1 2). For this particular case. ' 
the SRW becomes a fiat top profile where we can achieve nearly 100% light efficiency. For 
smaller and irregular phase patterns fine-tuning of tj in the region 0.4-0.6 provides for an 
efficient operating regime while maintaining minimal losses. 

In order to optimize a synthesized light distribution we wish to generate a dark background 
for the illumination at the observation plane as our principal design criterion, so we can write: 

^(*;.:vi;$o)=0 (20) 

where (jc^.j;^) indicates observation plane coordinates of the background of the generated 

illumination and ^0 is the relative phase shift generating a zero-Intensity level at the 
observation plane. 

Applying this dari< background condition In Eq. (13) we can obtain the following expression 
for a no-loss phase-only filter with filter transmission parameters. A = B = V. 

^|a|(l-e)qp(/e))=exp(4) (21) 

A key point arising from Eq. (21 ) is that we now have a simple way of expressing a new 
design criterion relating the spatial average value of any input phase pattem to the zero- 
order phase shift of a matched Fourier phase filter. Since K Is by definition positive and by 
taking the modulus of Eq. (21 ) we obtain: 

^|q = |2sin(e/2f (22) 

Eq. (22) is a key result for the fully transmissive wavefront engineered GPC mapping that 
makes it possible to deduce the range of valid phase parameters fulfilling our design criteria 
from Eq. (20). The largest possible value that the temi. K^ , takes on is unity, this leads to 
the following solution interval for Eq. (22) within a full phase-cycle: 

e = [n/3;57c/3] (23) 
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From Eq. (22) we also observe that k[a\ can take on a value limited to the Interval: 

^|a| = [X;l] (24) 

Eq. (22) and the solution Intervals described by Eqs. (23)-(24) specify the design parameters 
for achieving optimal performance In extracting and displaying the phase information earned 
by the incoming wavefront. Moreover. Eq. (22) hints towards extra means of optimisation by 
encoding the phase modulation depth itself In addition to the no-loss phase-only filter This 
last approach Is particularly useful when the filter phase has a restricted dynamic range or is 
fixed. 

Now. assuming that we have a fixed and fully transmissive phase-only filter, the best choice 
for the filter parameter is a value that allows for the largest dynamic range of phasor values 
at the input. Accordingly, the smallest possible /^a/ value. ira = 1/2. Is desirable Implying 
that e = K, leading to the output intensity distribution: 

/(y.y) = 2[l-cos(,i,(^',/))] (25) 

Inserting Ka^lll 1" Eq. (8) we obtain the following two requirements for the input 
encoded phase function 4{x,y) : 

KA-' IJcos (([) {x, y))dxdy = X 

KA-'\\sm(if{x,y))dxdy=.Q ^^^^ 

A 

We observe that It Is only the first requirement in Eq. (26) that is directly related to the output 
mtenslty in Eq. (25) via the cosine tem,. Since there are always two choices for a given 
phasor value that result in the same cosine value (excluding 0 and .). we notice that the 
second requirement can subsequently be fulfilled independently of the first requirement 
simply by complex conjugating an appropriate number of phasor values. This fact is a key 
feature of the GPC-method since it makes it possible to solely concentrate on the first 
requirement In the process of synthesizing a desired an virtually no-ioss grey level Intensity 
pattern. ^ 

The fl,^, requirement in Eq. (26) can be fuifiiied by several means, including: dynamic phase 
range adjustment, fill factor encoding, phase-hlstogram adjustment, spatial scaling of phasor 
pattern, raster encoding etc. In a histogram adjustment technique one will typically start out 
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wBh a desired naWve intensity distribution /(.'.y)^ where the maximum achievable 
intensity level is unknown but relative intensity levels a« known a^l the towest Intensity 
level ,s f«ed by the background criterion of Eq. (20). The pmcedure is now to adj^ the 
histogram for /(^.y)^ while maintaining IdenUca. relative Intensity level natfos until the 

by complex ooniugating an appropriate part of the phasors. The simples, p oLL 

ZZT ^ ^ "^'"^ "'^^^-*>y Of .he spatial 

^n. However, this -phasor fllppl^. procedure can also be turned Into an advantageous 

10 ^irr .he spatial frequency content In ordlto 

Ino^ 77T " ^' Plane by 

^ng «^e spat^i phasor location into account. E.g. neighbouring phasor values can be 

Chosen to have a maximum difference between them, thereby introducing high spatial 

frequency modulason easing the Altering in the spaUal Fourier domain. 

In most cases, however, equalized output Intensity levels are suflident. In the succeeding 

°" °' ^^-^ •« aohleverat 

«.t^ut .mens,ty ievete. A derivation based on ternary phase teveb allows for the widest 

iZlntr" ' ^"^ ^ »™ but 

important binary phase level encoding as a spedai case. 

For«,e ternary phase encoding, we considerthe i.umlnated portton of the input aperture 
20 -«a,Aasdi,.dedintosub^reas^^,«,d/.,withrespeotivephasevalues ^ and 
«>3 . We are aiming for the derivation of general expressions reiaSng the addressing 
parameterB for the phase modulation to the ,««e of possible phase parameters of Ihe 
Founer «ter obeying the des^n Cterion we have already set out We can express the total 
^ area and its average phase modutation. as the sum of «,e phase-w ighted ul 



(27) 



sucn tnat i^, = ^ / ^ and = ^2 / ^ : 

(l-i^-Fjexp(4) + irexp(4,) + 77exp(4) = |aj 
30 AS previously mentioned we are interested in binary Intensity patterns with levels 

corresponding to the Input phase values. In this case the dark background region Is defined 



(28) 
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by (4,.?o) and the bright output level of intensity. /. is determined by (4.^,) and J,) 
at the input plane. For the binary output intensity condition it follows that: 

^{*0 = 4*0 (29) 
This equality corresponds to a symmetric condition that can be easily verified applying the 
5 phasor chart analysis technique demonstrated In section 4. Due to this symmetry we can 
simplify the analysis by applying the following substitution: 

^^-h-h-k-k (30) 
so that Eq.(28) can be rewritten as: 

F.(exp(iA<|.)-l)+ir(exp(-fA<|,)-l) = i5:-(i_exp(/e))-'-l (g^j 

10 It is now a straightforward task to solve Eq. (31) for the real part and the imaginary part 
respectively, to obtain the following set of equations: 

+F, =(2i:-l)(2i:(l-cos(A<|))))"' 
F, -F, =sin(e)(2JS:sin(A(j,)(l-cos(e)))"' ^^^^ 



1 5 This can also be expressed in terms of the fractional 



areas, such that: 



^.=(4i^r[(2i:-l)(l-cos(A*))-Vsin(e)(sm(A<|,)(l-cos(e)))-'J 
^a=(4iS:r[(2ir-l)(l-cos(A4,))-'-sin(e)(sin(A<|>)(l-cos(e)^^^ 



(33) 



Since we have focused on solutions where identical intensity levels are obtained in both the 
Frregion and the F^reglon we can define the resulting illumination compression factor. C. In 
the following way: » ^. " » 

C.(/;+/r,)-=(,.(2js:)-)-(,_„,(^^)) (34) 

minimum compression factor con-asponds to uniform Illumination at tl,e output such that 
P,*F,-\. whereas the maximum compression factor is found to be C -> «. tor AT = X . 

25 An interesting special case can be deduced fremgq. (31) by setting J', - 0 . where we And 

that: 



wo 2004/113993 



-22- 



PCT/DK2004/000452 



F^F, =(^(l-exp(ie))-l)(i5:(l-exp(zA<|»))(l-exp(/e)))-' 05) 
implying that for the binary phase moduiation case we must have: 



~^ (36) 



in order for the fill factor, F. to be real-vaiued. 

This result turns out to be the special case that corresponds to the set of solutions where a 
binary phase pattern serves as the input. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 illustrates schematically a known phase contrast imaging system. 

Fig. 2 Illustrates schematically an embodiment of the present invention. 

Fig. 3 also schematically illustrates the embodiment of Fig. 2, 

Fig. 4 Illustrates various layouts of the spatial phase filter according to the present 
invention, and 

Fig. 5 Illustrates schematically another embodiment of the invention. 
DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
Fig. 1 shows a known phase contrast imaging system embodied in a 4f common path 
mterferometer 1. A laser emits a light beam, which is expanded by a beam expander into a 
Plane light wave of unifomr, intensity and directs it towards a phase modifying element The 
light beam is transmitted through the phase modifying element and a Fourier transforming 
lens (L1). The phase modifying element is positioned in the front focal plane of the lens (LI) 
and a spatial phase filter 6 is positioned in the back focal plane of the lens LI that is also the 
front focal plane of a lens L2. The Fourier transfomiing lenses LI, L2 need not have identical 
focal lengths. Different focal lengths lead to a magnification ratio different from one The 
spatial phase filter 6 phase shirts by 8, and optionally attenuates (by a factor B) the zero 
order diffraction part 8 of the light modulated by the phase modifying element. Optionally the 
remaining diffraction part of the light modulated by the phase modifying element may be' 
attenuated by a factor A. The reconstructed intensity pattern I(x'. y') Is generated In the back 
focal plane 9 of the lens L2. 

Figs. 2 and 3 illustrate schematically an embodiment 10 of the present invention in a 4f 
common path Interferometer configuration. An array 12 with a plurality of Individual light 

sources PS„ PS. PS„ emits a set of light beams Wave,. Wave. Wave„. which are 

colhmated by lens L3 into a plurality of plane light waves directed towards a phase modifying 
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element 14 to be Inddent upon the element at different respective angies of approach The 
ph^se modifying element 14 phase modulates the electromagnetic fields by phasor values 
e"^"'. (X. y) are coordinates of the resolution elements or pixels of the phase modifying 
element 14. The light waves are transmitted through the phase modifying element 14 and a 
Founer transforming lens LI. The phase modifying element 14 is positioned in the front focal 
plane of the lens L1 . and a spatial phase filter 6 is positioned in the back focal plane of the 
lens LI that is also the front focal plane of a lens L2. The Fourier transforming lenses LI L2 
need not have identical focal lengths. Different focal lengths lead to a magnification ratio* 
different from one. In each phase shifting region PCF,. PCF.. ... , PCF„. the spatial phase 

filter 6 phase shifts by respective values 9,. 9, 9, and optionally attenuates (by 

respective factors B,. B. B^) the respective n'th zero-order diffraction part of the light 

modulated by the phase modifying element 14. Optionally, the remaining diffraction part of 
the light modulated by the phase modifying element may be attenuated by a factor A The 
reconstructed Intensity pattern l(x'. y') Is generated In the back focal plane 9 of the lens L2. 
The illustrated system 10 may be controlled by a computer (not shown), comprising interface 
means for addressing each of the resolution elements of the phase modifying element 14 
and transmitting a phasor value e'^-v) to the addressed resolution element (x. y) Further 
the computer may comprise light control means for controlling the power of the llght-emitiing 
an-ay 12. The computer may also comprise input means, such as a keyboard, a diskette 
drive, an optical disc drive, a network interface, a modem, etc. for receiving an Image pattern 
to be synthesized by the system 10. From the received Image pattern, the computer may be 
adapted to calculate phasor values to be transmitted to the resolution elements (x. y) of the 
phase modifying element 14. e.g. based on a histogram technique as described herein 
Optionally, the phase shifts of the spatial phase filter 6 may be adjustable and controllable 
by optional phase control means of the computer, which may be further adapted to adjust 
the phase shift. 

Fig. 4 illustrates various layouts of the spatial phase filter according to the present invention 
for Illustration of the flexibility provided by the Invention. The illustrated layouts match 
corresponding layouts of source arrays. 

Fig. 5 Illustrates schematically an embodiment 10 of the present Invention in a 4f common 
path interferometer configuration. An array 12 with a layout 12a of a plurality of individual 
light sources 12b emits a set of light beams, which are collimated by lens L3 into a plurality 
of plane light waves directed towards a phase modifying element 14 to be incident upon the 
element 14 at different respective angles of approach. The phase modifying element 14 
phase modulates the electromagnetic fields by phasor values e'^^ v), (x. y) are coordinates of 
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the resolution elements or pixels of the phase modifying element 14. The light waves are 
transmitted through the phase modifying element 14 and a Fourier transforming lens LI The 
phase modifying element 14 is positioned in the front focal plane of the lens LI . and a spatial 
phase filter 6 with a layout 6a of a plurality of Individual phase shifting areas 6b is positioned 
in the back focal plane of the lens L1 that is also the front focal plane of a lens L2 The 
Founer transforming lenses L1 . L2 need not have identical focal lengths. Different focal 
lengths lead to a magnification ratio different from one. Each of the phase shifting areas 6b 
of the spatial phase filter 6 phase shifts light incident upon it by e„, and optionally attenuates 
(by a factor B). The light Incident upon it Is constituted by the zero order diffraction part 8 of 
light emitted from a conresponding laser 12b of the laser array 12a and modulated by the 
phase modifying element 14. Optionally, the remaining diffraction part of the light modulated 
by the phase modifying element may be attenuated by a factor A. The reconstructed 
intensity pattern Kx". y) is generated In the back focal plane 9 of the lens L2. 
The illustrated system 10 may be controlled by a computer 20 comprising Interface means 
for addressing each of the resolution elements of the phase modifying element 14 and 
transmitting a phasor value e'«-) to the addressed resolution element (x. y). Further the 
computer may comprise light control means for controlling the output of the an^y 12' The 
computer may also comprise input means, such as a keyboard, a mouse, a diskette drive 
an optical disc drive, a network interface, a modem, etc. for receiving an image pattern to be 
synthesized by the system 10. From the received image pattern, the computer may be 
adapted to calculate phasor values to be transmitted to the resolution elements (x. y) of the 
phase modifying element 14. e.g. based on a histogram technique as described herein 
Optionally, the phase shifts of the spatial phase filter 6 may be adjustable and controllable 
by optional phase control means of the computer, which may be further adapted to adjust 
the phase shift. An imaging system captures the synthesized intensity pattern (x" /) and 
transmits it to the computer 20 for display to the user and for possible automatic adjustment 
A user of the system may move individual light beams around using the mouse. 
PROPOSED APPLICATIONS: 

. Laser machine tool e.g. for shaping, forming, assembling, etc. such as cutting, drilling 
milling, planing, marking, branding, trimming, hardening, scribing, labelling, welding 
soldenng. two- and three-dimensional surfaces especially by use of CO2 and Nd YAG 
laser based systems. The main advantage is that energy is not absorbed in the system 
according to the present Invention (thereby preventing damage of the optical hardware) 
and this non-absorbed energy Is Instead utilized to Increase the intensity level of the 
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desired light distribution in the Image plane. High power can be delivered to selected 
regions on a work piece simultaneously. 

Optical tweezer array for manipulation of micro-objects, such as micro-components 
b-ological cells, etc. using electromagnetic gradient forces proportional to the optical 
intensity pointing in the direction of the intensity gradient. 

Efficient and dynamic spot-an-ay generators based on phase contrast Imaging. In order 
to provide bias or holding beams for anays of optoelectranic elements, such as bistable 
elements, photonic switches and smart pixels. 

Generation of structured light (loss less) for machine vision applications. E.g. periodic 
and skew periodic mesh grid illumination that can be updated In parallel. 
Photolithographic applications (laser 3D direct writing in parallel v>«thout the need for 
sequential scanning). E.g. high power laser direct writing of waveguides In Ge-doped 
silica. ^ 

Spatial light intensity modulation in general by use of pure phase modulation (radiation 

focusators). 

Laser beam shaping (dynamic). 

Highly efficient parallel image projection without the need for a laser-scanning device. 

Dynamic Infrared Scene Projection (DIRSP). 

Exposure device for grating and mask production. 

LIDAR applications. 

Laser printing in parallel. 

Laser show applications. 

Atmosphere research. 



